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We present a robust strategy to veto certain classes of instrumental glitches that appear at the 
output of interferometric gravitational-wave (GW) detectors. This veto method is 'physical' in 
the sense that, in order to veto a burst trigger, we make use of our knowledge of the coupling of 
different detector subsystems to the main detector output. The main idea behind this method is 
that the noise in an instrumental channel X can be transferred to the detector output (channel 
H) using the transfer function from X to H, provided the noise coupling is linear and the transfer 
function is unique. If a non-stationarity in channel H is causally related to one in channel X, the 
two have to be consistent with the transfer function. We formulate two methods for testing the 
consistency between the burst triggers in channel X and channel H. One method makes use of the 
null-stream constructed from channel H and the transferred channel X, and the second involves 
cross-correlating the two. We demonstrate the efficiency of the veto by 'injecting' instrumental 
glitches in the hardware of the GEO 600 detector. The veto safety is demonstrated by performing 
GW-like hardware injections. We also show an example application of this method using 5 days of 
data from the fifth science run of GEO 600. The method is found to have very high veto efficiency 
with a very low accidental veto rate. 

PACS numbers: 95.55.Ym, 04.80.Nn 
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I. INTRODUCTION 

The existence of gravitational waves (GWs) is the 
last and most intriguing prediction of Einstein's Gen- 
eral Theory of Relativity (GTR), which is yet to be 
verified by a direct observation. A worldwide net- 
work of GW detectors consisting of ground-based in- 
terferometers [l|, 0, [3, HI and resonant bars d, @, 0| , 
most of them operating at or close to their design 
sensitivity, has started looking for signatures of GWs 
produced by astrophysical and cosmological sources. 
Apart from providing excellent tests of GTR in the 
strongly-gravitating regions [1[, GW astronomy is ex- 
pected to open a new window to the universe providing 
unique information about various astrophysical phe- 
nomena. For example, GW observations from the final 
coalescence of neutron star binaries will shed light on 
the nuclear equation of state Q, and the GWs pro- 
duced by a neutron star or small black hole spiral- 
ing into a massive black hole will provide a 'map' of 
the spacetime geometry around the larger object |10| . 
Also, GW observations can be used to measure the 
Hubble constant-deceleration parameter, and cosmo- 
logical constant [II], [l2|, LL3| ■ Similarly, the observation 
of the stochastic GW background from the early uni- 
verse can test a number of current speculations about 
the very early universe [Til ]. 

The best understood (and, perhaps, the most 



'Electronic address: Ajith.Paramcswaran@aei.mpg.dc 



promising) astrophysical sources of GWs for ground- 
based interferometers are the inspiralling compact bi- 
naries consisting of black holes and/or neutron stars. 
What makes this class of sources extremely interesting 
is that the expected waveforms can be very accurately 
modelled (and easily parametrized by a few intrinsic 
parameters of the binary, like the component masses 
and spins) by approximation techniques in GTR. This 
allows the data analysts to use the matched filtering 
technique in order to extract the signal that is buried 
in the detector noise. Another important class of astro- 
physical sources which can be probed through optimal 
filtering technique is spinning neutron stars and pul- 
sars. 

There are other classes of GW sources, like core- 
collapse of massive stars in supernovae, binary black 
hole/neutron star mergers, accretion induced collapse 
of white dwarf stars, Gamma ray bursts etc., for which 
the physics is largely unknown, or too complex to al- 
low computation of detailed gravitational waveforms. 
These kind of sources are generally classified as 'un- 
modelled burst sources'. GW observations from merg- 
ing black hole binaries can potentially bring very use- 
ful insights to the nonlinear dynamics of the spacetime 
curvature as the two black holes convert themselves 
to form a single black hole. GWs from a Gamma ray 
burst are expected to carry detailed information about 
its source, which can not be probed via electromag- 
netic observations. Similarly, correlated neutrino and 
GW observations from a core-collapse supernova could 
bring interesting insights into the newborn neutron 
star/black hole in the core of the supernova (See [Hj] 
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for a detailed review). 

Any possible GW signal coming from an astrophys- 
ical/cosmological source is generally buried in the de- 
tector noise. To extract these GW signatures from the 
noise is a nontrivial data analysis problem. Indeed, if 
time evolution of the GW phase is accurately known, 
the optimal filter for searching for this signal buried in 
the noise is the well known matched-filter. But, since 
matched filtering relies on the prior knowledge of the 
signal, it may not be the best detection strategy in the 
search for unmodelled, short-lived GW bursts. In the 
next section, we give a brief introduction to the data 
analysis for 'unmodelled burst sources'. 



A. The search for transient, unmodelled 
gravitational-wave bursts 

One class of search methods that is being employed 
in the burst data analysis is based on time-frequency 
decomposition of detector data. These algorithms con- 
struct time-frequency maps of the time-series data 
and look for 'time- frequency regions' containing ex- 
cess power which are statistically unlikely to be as- 
sociated with the background noise distribution (IH ]. 
Some of these algorithms are based on clustering the 
'time-frequency pixels' containing excess power, and 
applying another threshold on these clusters of pix- 
els [17] . Time- frequency detection algorithms using 
basis functions other than the standard Fourier ba- 
sis functions are also proposed [H, EH ■ Another class 
of burst detection algorithms look for slopes or ridges 
in the time-series data, or in its time-frequency repre- 
sentation [lO, HH, HH. In general, these methods are 
claimed to be robust in detecting short-lived signals 
with minimum a priori information. 

While the optimal filtering technique, along with ac- 
curate models of the waveforms, enables one to accu- 
rately estimate the physical parameters (such as masses 
and spins) of the GW source, the time-frequency meth- 
ods, by construction, are unable to accomplish this. 
Instead, these algorithms try to parametrize the un- 
derlying gravitational waveforms using a set of quan- 
tities like the characteristic central frequency, dura- 
tion, bandwidth etc. The detection algorithms, im- 
plemented in the data analysis pipelines, are usually 
referred to as event trigger generators (ETGs). 

Since current interferometric detectors are highly 
complex instruments, the detector output typically 
contains a large number of noise transients, or 
'glitches', of instrumental origin which cause the ETGs 
to generate spurious triggers. One of the main chal- 
lenges in the burst data analysis is to distinguish these 
spurious bursts from actual GW bursts. Since the ex- 
pected GW signals are unmodelled, it is practically 
impossible to distinguish these 'false alarms' from ac- 
tual GW bursts based on their signal characteristics. 
One way of dealing with this issue is to require that the 
triggers be coincident (within a time window) in mul- 
tiple detectors located at different parts of the world. 



Although this 'coincidence requirement' reduces the 
list of candidate triggers by a considerable amount, 
this does not completely cure the problem. While 
coincident instrumental bursts in multiple detectors 
are highly improbable, long data-taking runs (typi- 
cally several months long) using multiple detectors can 
produce a large number of random coincidences (po- 
tentially thousands per month [HI). It is thus very 
important to develop robust techniques to distinguish 
between true GW bursts and spurious instrumental 
bursts which are coincident in different detectors - pop- 
ularly known as veto techniques. 

Since a number of environmental and instrumental 
noise sources can potentially couple to the main detec- 
tor output, many such noise sources are continuously 
recorded along with the data from the main detector 
output. The measurement points for time-series data 
within the detector are referred to as 'channels'. One 
class of veto methods is based on identifying triggers 
in the 'gravitational-wave channel' (the main detector 
output) which are coincident with triggers in an in- 
strumental/environmental noise channel. The 'coinci- 
dence windows' are chosen such that the 'accidental' 
(random) coincidence rate between the two channels is 
limited to an acceptable amount. See [23|,[24],[25|,[26| for 
some recent work on such 'statistical vetoes'. Another 
class of 'physical vetoes' is based on our understanding 
of how a GW should (or, should not) appear in certain 
channels [27], [H, [29[ . Moreover, a number of 'waveform 
consistency tests' between multiple detectors are also 
employed in the burst searches [3(3, HH, [H, HH, [H| ■ 

In this paper, we demonstrate a veto strategy which 
makes use of our understanding of the physical cou- 
pling of various detector subsystems to the detector 
output. This method is different from the physical veto 
methods discussed above in the sense that, here we use 
the knowledge of the coupling mechanism involved in 
transporting glitches from a particular subsystem to 
the main detector output. The main idea behind this 
method is that the noise in an instrumental channel 
X can be transferred into the detector output (chan- 
nel H) using the transfer function from X to H 1 . If 
a noise transient in channel H is causally related to 
one in channel X, the two have to be consistent with 
the transfer function. The basic idea of this method 
was formulated in [35], which also demonstrated this 
in the presence of stationary Gaussian noise. Here, we 
demonstrate a more general formulation which admits 
non-Gaussian tails in the noise distribution and other 
'real- life effects' (see Sec. Ill Cp . We also propose an 
alternative statistic to test the consistency of the two 
triggers (see Sec. Ill B[) . Sec. HI1 reviews the basic ideas 
of this method. In Sec. lIIII we demonstrate the method 



Throughout this paper, channel X refers to the measure- 
ment point for time-series data from a detector subsys- 
tem/environmental noise source X, and channel H refers to 
the main detector output (the 'GW channel'). 
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by performing hardware injections in the GEO 600 de- 
tector 0] . The results of the veto analysis performed on 
5 days of data from the fifth science run of GEO 600 
are discussed in Sec. IIVI In Sec. [Vl wc summarise our 
main findings. 



II. VETOES USING KNOWN 
INSTRUMENTAL COUPLINGS 

Let Xi and hi denote the discretely sampled time- 
series data recorded in channel X (which records the 
noise from a detector subsystem) and channel H (the 
'gravitational-wave channel'), respectively. We denote 
the corresponding discrete Fourier transforms (DFT) 
by Xk and hk, respectively. The input and output of a 
linear, time-invariant system are related by the transfer 
function Tk of the system, which is defined as 



pXH 
jtXH _ r k 
k ~ pXX ' 



(2.1) 



Here P fc XH - ~ 



Xkht is the cross-power spectral density 
of Xi and hi, and -P* x = Xkxl is the power spectral 
density of Xi, where the 'bars' indicate ensemble aver- 
ages. If the coupling of noise between channel X and 
H is linear and the transfer function is time- invariant, 
the Fourier transform of the noise measured in chan- 
nel X at any time can be transferred to channel H, by 
using the transfer function 



Xk T, 



XH 



(2.2) 



Xk and hk can be thought of as components of two 
vectors x and h defined in two TV-dimensional linear 
vector spaces. In the mathematical sense, Eq. (|2.2p 
maps x to the space of h. In the physical sense, this 
means that if a noise transient originates in X , one can 
predict how it will appear in H. If there exists a noise 
transient in H at the same time 2 , we can compare it 
with the above 'prediction'. If a noise transient in the 
channel H is causally related to one in channel X, the 
data vectors x and h have to be consistent with the 
transfer function. This allows us to formulate a power- 
ful strategy to veto noise transients originating within 
the detector. 

The basic idea is the following: firstly, we identify 
time-coincident burst triggers in channels X and H. 
We compute the DFTs of two short segments of data 
in channels X and H. The length of these segments 
(typically a few tens of milliseconds) is chosen so as 
to encompass only the noise transient under investiga- 
tion. If these two noise vectors are consistent with the 




FIG. 1: Schematic diagram of the linear vector space in 
which the analysis methods are formulated. 



transfer function, as given by Eq. (|2.ip , it is highly likely 
that the noise transient originates in X, and we veto 
the trigger. On the other hand, if the two noise vec- 
tors are inconsistent with the transfer function, we con- 
clude that this particular noise transient in H does not 
have its origin in X. In the following subsections, we 
construct two different statistics which can be used to 
make this decision. One statistic is based on construct- 
ing a null-stream between channel H and the 'trans- 
ferred' channel X; i.e., if the noise transient originates 
in subsystem X, and is sensed by channel X, it is pos- 
sible to construct a linear combination of the two data 
streams such that it does not contain any excess power. 
The second statistic is based on the cross-correlation 
of the noise in channel H with the 'transferred' noise 
in channel X. 

In the following two subsections, we assume that 
the data streams Xi and hi are drawn from zero- 
mean Gaussian distributions. Also, we assume that 
the transfer function from X to H is accurately mea- 
sured/calculated and is time-invariant. In the third 
subsection, we discuss the 'real-life scenario' where the 
data streams are not perfectly Gaussian and the trans- 
fer function is non-stationary. 

The linear vector space in which the analysis meth- 
ods arc formulated is schematically illustrated in Fig.[T] 



A. Using the null-stream 

The null-stream between h and x' is the compo- 
nent of the vector h orthogonal to x'. This can be 
constructed using the Gram-Schmidt orthogonalisa- 
tion [H: 

<5 = h - proj x ,h, (2.3) 
where we define the projection operator by 



proj a v 



u, u 



(2.4) 



2 The time-coincidence window should be chosen according to 
the typical time scale of the transients that we are concerned 
with. We use a time window of a few tens of milliseconds 
since the current searches for GW bursts seek to detect bursts 
of duration << 1 sec. 



In the above expression (v, u) denotes the inner prod- 
uct between the vectors v and u: 



v,u 



v k u* k . 



(2.5) 
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If the noise transient originates in X, it will completely 
disappear in d. In order to test this, we compute the 
excess power statistic [Til ] from 5: 



es 



E 



\5 k 



(2.6) 



where o~\ is the expected variance of 5k in the absence of 
any excess power. In the absence of any excess power in 
8, es will follow a Gamma distribution [3{|. The scale 
parameter a and shape parameter (3 can be estimated 
from the stationary noise (see [35[ for more details). 

If es is less than, or equal to, a threshold r, we veto 
the trigger. The threshold r giving a rejection proba- 
bility <f> (probability that a 'causal' trigger is vetoed) 
can be calculated from 



<I> 



T(x; a, f3) dec, 



(2.7) 



where T(x]a,f3) is the probability density of the 
Gamma distribution with parameters a and (3. 
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FIG. 2: Schematic illustration of the veto pipeline. 



B. Using the cross-correlation 



C. Implementation 



The linear cross-correlation coefficient between two 
vectors x' and h is the cosine of the angle between 
them: 



Re- 



x',h 



(2.8) 



where ||u|| denotes magnitude of the vector u. If the 
noise transient in channel H indeed originates in X, 
x' and h should display a high correlation. On the 
other hand, if the noise transient does not originate in 
X , the vector x' and h will be randomly oriented, and 
hence the linear cross-correlation coefficient r will tend 
to be small in absolute value. This can be converted 
to the normally distributed variable z by the Fisher 
transformation 1371: 



1 



■In 



1 



1 -r 



(2.9) 



The new variable z will be approximately normally 
distributed with mean zero and standard deviation 
\j\JX — 3, where N is the dimension of the vectors 
x' and h. 

If z is greater than, or equal to, a threshold A, we 
veto the trigger. The threshold giving an accidental 
veto probability of ip can be calculated from 



f(x;n,a 2 )dx, 



(2.10) 



The assumption we made in the previous subsec- 
tions that the transfer function is time invariant is 
strictly not true. Transfer functions in actual detec- 
tors can vary in time. The slow temporal variation 
of the transfer function can be taken into account by 
making repeated measurements of the transfer func- 
tion and tracking its evolution by continuously inject- 
ing and measuring spectral lines at certain frequencies 
(sec [38]). But the non-stationarities of the transfer 
function on short time scales are hard to track. It may 
also be noted that the noise in the present-generation 
interferometers is not stationary Gaussian and exhibits 
tails in the distribution. Considering these 'real-life' 
effects, it may not be wise to use the 'ideal-case' rela- 
tions given by Eqs. (|2.7[ ) and (|2.10[) to compute the veto 
thresholds. For instance, due to the imperfect transfer 
function, the 'null-stream' 6 can contain some 'residual 
burst', and, as a result, the excess power statistic com- 
puted from S will not fall into the expected Gamma dis- 
tribution. But, we do expect the excess power statistic 
es computed from 8 to be smaller than the same (eh) 
computed from h. If the ratio s = e^/es is greater than 
a threshold, we veto the trigger. The veto threshold 
corresponding to a certain accidental veto probability 
is calculated as described below. 

We time shift Xi with respect to hi to destroy the 
causal relationship between the two data streams 3 . 
The coincident triggers in the time-shifted data 



where f(x;/j,,cr ) is the probability density of the nor- 
mal distribution with mean /i = and variance a 2 = 
l/(JV-3). 



3 Time shift analysis is commonly employed in burst searches in 
order to estimate the accidental consistency, or 'background' 
rate. See, for example, [23l |. 
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streams arc identified and the 'excess-power ratio', s, 
for each coincident trigger is calculated, n such time 
shifts are performed to get better statistics from the 
data. A threshold, t s , is chosen such that only an 
acceptable number of coincident triggers in the time- 
shifted analysis have s > r s . This threshold r s is used 
to veto the triggers in the 'zero-lag' analysis (without 
time shifting the data). The time-shifted analysis can 
also be used to calculate the veto threshold A for the 
analysis using the cross-correlation statistic. Here, we 
choose a threshold A such that only an acceptable num- 
ber of coincident triggers have z > A in the time-shifted 
analysis and use this threshold to do the zero-lag anal- 
ysis. The veto pipeline is schematically illustrated in 
Fig.H 



D. Caveats 

It is worth stressing that this method relies on the 
linearity in the coupling of the noise from the detec- 
tor subsystem X to the detector output, and can not 
be used where the coupling is nonlinear. This method 
also assumes that the transfer function from X to H 
is unique, and channel X accurately senses the distur- 
bances in X. In other words, this technique can be 
only applied to systems that exhibit a linear coupling 
through a set path, or multiple paths that are fixed. 
An environmental monitor will often fail to meet this 
requirement, unless the sensors are exceptionally well 
placed, because each local disturbance could couple dif- 
ferently into the monitor and the GW channel, mean- 
ing that a different transfer function would be needed 
for each physical point of origin for the disturbance. 



III. ANALYSIS ON HARDWARE-INJECTED 
BURST SIGNALS 

A. Injections mimicking instrumental bursts 

In order to test this veto method, we injected around 
300 sine-Gaussian burst signals over a period of one 
hour into four subsystems of GEO 600, whose cou- 
plings to H were known and well understood; the in- 
jections were performed serially, one subsystem after 
another. The four subsystems we chose are listed be- 
low. These descriptions are technical and concise. For 
more information, refer to [io| . 

Laser amplitude noise (P rc f): We make bursts of laser 
amplitude noise by injecting glitches into the laser am- 
plitude stabilisation loop. We detect these glitches by 
measuring the light power reflected from the Power- 
recycling cavity in the data acquisition system as chan- 
nel P re f. 

Laser frequency noise (F TC {): We make bursts of laser 
frequency noise by adding glitches to the error-point 
of the Michelson common-mode control servo which 
keeps the Power-recycling resonant by adjusting the 
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FIG. 3: Time shift analysis on 'instrumental-glitch-like' 
hardware injections performed in the $mi channel. The 
horizontal axis shows the time shift applied between Xi and 
hi, and the vertical axis shows the excess-power ratio s. 
The black dots correspond to the coincident triggers in the 
time-shifted analysis and the red crosses correspond to the 
ones in the zero-lag analysis. 
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FIG. 4: Same as in Fig. [3] except that the vertical axis 
shows the cross-correlation statistic z. The black dots cor- 
respond to the coincident triggers in the time-shifted anal- 
ysis and the blue crosses correspond to the ones in the zero- 
lag analysis. 



frequency of the master laser. The recording of this 
error-point in the data acquisition system serves as the 
veto channel, F re {. 

Michelson oscillator phase noise ($mi): The Michelson 
differential arm-length in GEO 600 is controlled by im- 
posing phase-modulation side-bands on the light enter- 
ing the interferometer. A heterodyne readout scheme 
is then used to derive an error signal which is fed back 
to the end mirrors of the Michelson to keep it on a dark 
fringe. We make glitches in the phase of the oscillator 
used to add the modulation sidebands by driving the 
voltage-frequency-control input of the crystal oscillator 
used to create this modulation signal. We phase-lock 
a reference crystal oscillator to the main crystal os- 
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cillator and the error-point of the phase-locked loop, 
which is sensitive to phase fluctuations on both oscilla- 
tor signals, is recorded in the data acquisition system 
as $mi and serves as a sensitive measurement of the 
phase noise on the main Michelson modulation side- 
bands. 

Michelson oscillator amplitude noise (Ami)' The am- 
plitude of the the main crystal oscillator is also sta- 
bilised to a quiet DC reference. We can add signals to 
the error-point of this stabilisation servo so as to im- 
pose additional amplitude noise on the main Michelson 
modulation signal. We added glitch signals in to this 
control loop and recorded its error-point in the data 
acquisition system as Ami to serve as a veto channel. 

Burst triggers in the veto channel and the GW chan- 
nel are generated using the mHACR [25|, HH burst 
detection algorithm. mHACR belongs to the class 
of time-frequency detection algorithms those make a 
time-frequency map of the data and identify time- 
frequency pixels containing excess power which are sta- 
tistically unlikely to be associated with the underlying 
noise distribution. For a detailed description of the al- 
gorithm and its performance, see (25[. However, we re- 
mind the reader that the details of the burst detection 
algorithm are immaterial as far as this veto method is 
concerned. The burst ETG is only used to identify the 
coincident triggers in the two channels, and any ETG 
with proper time estimation of the burst event should 
serve this purpose. 

Coincident triggers in the two channels arc identi- 
fied using a time window of ±10 ms. The results of 
the time shift analysis on the hardware injections per- 
formed in the <I>mi are shown in Figs. [3] and [4] The 
horizontal axis shows the time shift applied between Xi 
and hi. The vertical axis in Fig. [3] shows the excess- 
power ratio s = e^/eg. The (black) dots correspond to 
the coincident triggers in the time-shifted analysis and 
the (red) crosses correspond to the ones in the zero-lag 
analysis. From this, a veto threshold of 2.35 is cho- 
sen which corresponds to an accidental veto rate of 1 
per day. All the coincident triggers in the zero lag are 
vetoed using this threshold. The vertical axis in Fig. 
S] shows the (transformed) cross-correlation statistic z. 
The (black) dots correspond to the coincident triggers 
in the time-shifted analysis and the (blue) crosses cor- 
respond to the ones in the zero-lag analysis. The veto 
threshold corresponding to an accidental veto rate of 
1 per day is 0.27, which resulted in vetoing 99% of the 
coincident triggers in the zero-lag. 

The veto analysis is performed on all the four chan- 
nels in which the hardware injections are done. Results 
of the analysis are summarised in Table |TJ It can be 
seen that only ~ 5% of the coincident triggers in the 
time-shifted analysis are vetoed, while more than 90% 
of the coincident triggers in the zero-lag are vetoed. 
This implies that the accidental rate of the veto is only 
~ 5% of that of the standard statistical veto (using 
a time window of ±10 ms) for almost the same veto 
efficiencies. 
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-frcf 


2.51 0.33 


4.48 x 10" 2 


0.90 


0.90 


Prof 


1.94 0.23 


5.45 x 10~ 2 


1.00 


1.00 




2.35 0.27 


6.12 x 10~ 2 


1.00 


0.99 


^Imi 


1.50 0.26 


4.62 x 10" 2 


0.97 


0.97 



TABLE I: Summary of the veto analysis on hardware in- 
jections mimicking instrumental bursts. r s and A are the 
chosen veto thresholds on the excess-power ratio s and the 
cross-correlation statistic z, respectively. £ is the fraction of 
coincident events that are vetoed in the time-shifted analy- 
sis. The fraction of coincident events vetoed in the zero-lag 
using the s statistic is denoted by Xs, while the same us- 
ing the 2 statistic is denoted by Xz • The chosen thresholds 
correspond to an accidental veto rate of 1 per day. 



B. Injections mimicking gravitational-wave 
bursts 

Some of the interferometer channels are sensitive to 
GWs to some non- negligible level. This raises the ques- 
tion of veto safety while using interferometer channels 
as veto channels, i.e., we have to make sure that we do 
not veto actual GW bursts which are coincident in the 
two channels. We argue that, since actual GW bursts 
are not causally related to the instrumental channels, 
the coincident triggers in channels X and H will not be 
consistent with the transfer function from X to H and 
hence, will not be vetoed using this method. Although 
the four channels under investigation in this paper are 
not expected to show any non-negligible sensitivity to 
GWs, there can be unexpected couplings, for example, 
through electrical faults or cross couplings in the data 
acquisition system. It is therefore prudent to explicitly 
demonstrate the safety of this veto method by doing 
GW-likc hardware injections. 

Hardware injections are performed by injecting sig- 
nals into the electrostatic actuators used to control the 
diffcrcntial-arm-lcngth degree of freedom of GEO 600. 
For the test described here, around 300 sine-Gaussian 
bursts were injected with varying amplitudes and with 
central frequencies in the range 200 to 1300 Hz. 

Figs. O and [5] show the results of the veto analysis 
performed on the GW-likc hardware injections. _F re f is 
used as the veto channel. It can be seen that neither of 
the test statistics (s or z) in the zero-lag analysis shows 
any excess significance over the corresponding time- 
shifted analysis. The veto thresholds corresponding 
to an accidental rate of 1 per day are r s = 1.64 and 
A = 0.35. Using these thresholds, we don't veto any of 
the injections. 



IV. AN EXAMPLE APPLICATION 

GEO 600 participated full time in the 'fifth science 
run' (S5 run) in coincidence with the LIGO detectors 
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FIG. 5: Time shift analysis on the GW-like hardware injec- 
tions. jf? ro f is used as the veto channel. The horizontal axis 
shows the time shift applied between xt and hi, and the ver- 
tical axis shows the excess-power ratio s. The black dots 
correspond to the coincident triggers in the time-shifted 
analysis and the red crosses correspond to the ones in the 
zero-lag analysis. 



FIG. 7: Time-shifted analysis on 5 days of data from the 
Fifth science run of GEO 600 using F le f as the veto chan- 
nel. The horizontal axis shows the time shift applied be- 
tween Xi and hi, and the vertical axis shows the excess- 
power ratio s. The black dots correspond to the coincident 
triggers in the time-shifted analysis and the red crosses cor- 
respond to the ones in the zero-lag analysis. 
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FIG. 6: Same as in Fig. [5] except that the vertical axis 
shows the cross-correlation statistic z. The black dots cor- 
respond to the coincident triggers in the time-shifted anal- 
ysis and the blue crosses correspond to the ones in the zero- 
lag analysis. 
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FIG. 8: Same as in Fig. [7] except that the vertical axis 
shows the cross-correlation statistic z. The black dots cor- 
respond to the coincident triggers in the time-shifted anal- 
ysis and the blue crosses correspond to the ones in the zero- 
lag analysis. 



from May 2006 to October 2006. The first few weeks of 
H data contained an additional population of glitches 
identified as coming from the laser frequency stabil- 
isation control loop. These excess glitches had cen- 
tral frequencies typically around 2 kHz. In fact, the 
glitches were broad-band in the frequency stabilisation 
loop, but the coupling of frequency noise to H is most 
prominent around 2 kHz and so this is where we see the 
excess noise in H. The identification and repair of the 
source of these glitches took several weeks. So we iden- 
tified an appropriate measure (-F re f) of the frequency- 
noise glitches that could be used as a veto channel (see 
Sec. IIII Al to see how this channel is derived). 

Veto analysis is performed on 5 days of data from the 
period described above (in the frequency range of 400 



Hz - 2kHz). Burst triggers in the two channels are gen- 
erated by the mHACR burst-detection algorithm. Co- 
incident triggers within the two channels are identified 
using a time window of ±10 ms for time-coincidence. 
Only triggers with signal-to-noise ratio > 6 are con- 
sidered for this analysis. Out of 5326 triggers in the 
GW channel, 2048 triggers were found to be coinci- 
dent with the i^ ro f channel. The accidental rate of 
the veto is estimated by doing 76 time shifts (from 
-100s to 100s). Fig. [7] shows the excess-power ratio s 
computed from the coincident triggers plotted against 
the applied time shift between the data streams, while 
Fig. [5] shows the cross-correlation statistic z plotted 
against the time shift, the dots correspond to the co- 
incident triggers from the time-shifted analysis and the 
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FIG. 9: Histograms of the excess-power ratio s computed 
from the time-shifted analysis (left) and the zero-lag anal- 
ysis (right). 
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FIG. 10: Histograms of the cross-correlation statistic z 
computed from the time-shifted analysis (left) and the zero- 
lag analysis (right). 



crosses correspond to those from the zero-lag analysis. 
We choose an accidental veto rate of 1 per day. The 
thresholds on the two statistics are estimated from the 
time-shifted analysis. This corresponds to a threshold 
of t s = 2.25 for the excess-power ratio s and a threshold 
of A = 0.54 for the cross-correlation statistic z. In the 
analysis using the null-stream, all coincident triggers 
with s > t s are vetoed, while in the analysis using the 
cross-correlation statistic, all triggers with z > A are 
vetoed. It was found that 88% of the coincident trig- 
gers are vetoed using the null-stream method and 92% 
of the coincident triggers are vetoed using the cross- 
correlation method. These correspond to 34% and 35% 
of the total number of H triggers in the data. 

Histograms of the two test statistics s and z com- 
puted from the coincident triggers are plotted in Figs. [5] 
and [TO] The plots on the left show the distributions of 
the test statistics computed from the time-shifted anal- 
ysis, normalised by the number of time shifts applied. 
These are the expected distributions of s and z in the 
absence of any causal relation between triggers in X 
and H (for the given data set) . Histograms on the right 
show the distributions of s and z computed from the 
zero-lag analysis. Two different populations are clearly 
visible in these plots. One population (centered around 
1 in the histograms of s; centered around zero in the 
histograms of z) corresponds to the triggers which are 
accidentally time-coincident in the channels X and H , 
while the other population (centered around 6 in the 
histogram of s; centered around 1.6 in the histogram of 



FIG. 11: A time-frequency plot of mHACR triggers from 
5 days of GEO 600 data. The horizontal axis shows the 
time and the vertical axis shows the frequency of the the 
burst triggers in channel H, as estimated by mHACR. Only 
those triggers which are coincident with F IC { are plotted. 
The (green) circles correspond to coincident triggers which 
are vetoed using the null-stream method, and the (black) 
dots correspond to the ones which are vetoed using the 
cross-correlation method. Coincident triggers which are not 
vetoed by any of the methods are indicated by (red) crosses. 
The chosen veto thresholds correspond to one accidental 
veto per day. 



z) corresponds to triggers in H which are causally re- 
lated to the ones in A". It is interesting to note that this 
method is able to distinguish clearly between these two 
populations. The reader may note that the number of 
accidental coincidences in the time-shifted analysis is 
~ 35% larger than that in the zero-lag analysis. This 
can be explained in the following way: the veto analysis 
has shown that around 35% of the triggers in channel 
H are causally related to F Te {. These 'causal' triggers 
will fall into the populations on the right in the zero- 
lag analysis (centered around 6 in the histograms of s; 
centered around 1.6 in the histograms of z). But, since 
the accidental coincidence rate is directly proportional 
to the total number of triggers, the presence of these 
causal triggers in the data would increase the coinci- 
dence rate in the time-shifted analysis by ~ 35%, thus 
explaining the excess-coincidences that we observe. 

Fig. QT] shows a time-frequency plot of the mHACR 
triggers from 5 days of GEO 600 data. The green cir- 
cles correspond to the coincident triggers (in channel 
H) which are vetoed using the null-stream method and 
the black dots correspond to the ones which are vetoed 
using the cross-correlation method. The red crosses 
correspond to the coincident triggers which are not ve- 
toed by any of the methods. 

A summary of the analyses performed using differ- 
ent accidental veto rates are given in Table [TT1 Also, in 
Fig. [T21 we plot the fraction of coincident events that 
are vetoed in the zero-lag (a measure of the efficiency 
of the veto) against the fraction of the coincident trig- 
gers which are vetoed in the time-shifted analysis (a 
measure of the accidental veto probability). This plot 
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Fraction of coincident events vetoed in time-shift 

FIG. 12: Fraction of coincident triggers that are vetoed in 
the zero-lag plotted against the fraction of the coincident 
triggers that are vetoed in the time-shifted analysis. The 
solid curve corresponds to the analysis using null-stream 
and the dashed curve corresponds to the analysis using 
cross-correlation. The triangles, the squares, and the dots 
correspond to accidental veto rates of 1 per month, 1 per 
week, and 1 per day, respectively. 



can be thought of as a receiver operating characteris- 
tic [il[ plot for this analysis, and can be used to choose 
thresholds which correspond to acceptable values of 
veto efficiency and accidental veto rate/probability. In 
the Figure, the solid curve corresponds to the analy- 
sis using null-stream and the dashed curve corresponds 
to the analysis using cross-correlation. It can be seen 
that, for high values of accidental veto probability, the 
two methods perform equally well. But for low values 
of accidental veto probability (< 2 x 10 -2 ), the curve 
corresponding to the null-stream analysis starts to fall 
off, and the cross-correlation analysis continues to per- 
form well. 



Accidental 


Threshold 


Veto fraction 


rate 


T s X 


£ Xs 




day" 1 


2.25 0.54 


1.73 x 10~ 2 0.88 


0.92 


week - 


5.09 0.86 


2.47 x 10~ 3 0.45 


0.90 


month -1 


7.40 1.11 


5.93 x 10" 4 0.23 


0.85 



TABLE II: Summary of the veto analysis on 5 days of data 
from GEO 600. r s and A are the chosen veto thresholds on 
the excess-power ratio s and the cross-correlation statistic 
z, respectively. £ is the fraction of coincident events that 
are vetoed in the time-shifted analysis. The fraction of 
coincident events vetoed in the zero-lag using the s statistic 
is denoted by Xs, while the same using the z statistic is 
denoted by Xz- The chosen veto thresholds correspond to 
the accidental veto rates tabulated in the first column. 



Fig. \T3\ provides a rough comparison between the 
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FIG. 13: Excess-power ratio s (horizontal axis) computed 
from the coincident triggers plotted against the cross- 
correlation statistic z (vertical axis). The vertical (red) and 
horizontal (blue) lines correspond to the veto thresholds t s 
and A on the two statistics, respectively. 

abilities of the two test statistics (s and z) in vetoing 
the instrumental glitches. The horizontal axis shows 
the excess-power ratio s and the vertical axis shows the 
cross-correlation statistic z computed from the coinci- 
dent triggers. The vertical (red) and horizontal (blue) 
lines in the plot correspond to the veto thresholds t s 
and A on the two statistics, respectively. Triggers on 
the right of the vertical line are vetoed by s, and those 
above the horizontal line are vetoed by z. 33.5% of the 
total number of H triggers are vetoed by both meth- 
ods. There exists a small population (~ 1% of the total 
number of H triggers) which is vetoed by z; but not 
by s, which suggests that z is a more sensitive statistic 
than s. But this may not be taken as a general in- 
dication that the cross-correlation is a more sensitive 
method than the null-stream. One can construct al- 
ternative statistics using the null-stream, which could 
be more sensitive than s. One possible alternative is 
e<5/eopt, where e opt is the excess power statistic com- 
puted from the optimal combination [42|, EH of the 
noise vectors h and x'. 



V. SUMMARY AND OUTLOOK 

One of the most challenging problems in the search 
for unmodelled GW bursts using ground-based detec- 
tors is to distinguish between actual GW bursts and 
spurious instrumental bursts that trigger the detection 
algorithms. In this paper, we proposed a veto method 
which makes use of the information on the physical cou- 
pling of different detector subsystems to the main de- 
tector output. We also demonstrated this method us- 
ing the data of the GEO 600 detector. By performing 
hardware injections mimicking instrumental glitches, 
we showed that glitches originating in a detector sub- 
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system can be vetoed using the transfer function from 
the subsystem to the detector output. We also ad- 
dressed the issue of veto safety by performing hard- 
ware injections mimicking GW bursts into GEO 600, 
and by showing that such injections are not vetoed. Fi- 
nally, we used this strategy to veto glitches in the data 
from the fifth science run of GEO 600, using the laser 
frequency noise channel as the veto channel. The anal- 
ysis was performed on 5 days of GEO 600 data from 
the second month of the science run. 35% of the 5326 
triggers in the GW channel were vetoed with an ac- 
cidental rate of 1 per day using the 'cross-correlation' 
method, while 34% of the triggers were vetoed using 
the 'null-stream' method. 

The method relies on linearity in the coupling of 
the noise from a detector subsystem to the detec- 
tor output, and the measurability/calculability and 
uniqueness of the transfer function. The assumption 
of linear coupling is valid as far as many detector 
subsystems in the large-scale interferometers are con- 
cerned. Strictly speaking, this method also requires 
time-invariant transfer functions. The way to track 
down slow temporal variations in the transfer func- 
tions is discussed in the paper. The formulation that 
we have developed was found to be robust against non- 



stationarities of short time scales in the transfer func- 
tions, and non-Gaussian tails in the noise distribution. 

When possible, using physical information has clear 
advantages over relying only on statistical correlations. 
The method that was proposed in this paper is a fully 
coherent way of testing the consistency of the glitches 
in the GW channel with those in an instrumental chan- 
nel. The authors hope that this will serve as a first step 
for developing a class of 'physical instrumental vetoes' 
for present and future detectors. 
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